Wehrl entropy, Lieb conjecture and entanglement monotones 
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We propose to quantify the entanglement of pure states oi N x N bipartite quantum systems 
by defining its Husimi distribution with respect to SU{N) x SU{N) coherent states. The Wehrl 
entropy is minimal if and only if the analyzed pure state is separable. The excess of the Wehrl 
entropy is shown to be equal to the subentropy of the mixed state obtained by partial trace of the 
bipartite pure state. This quantity, as well as the generalized (Renyi) subentropies, are proved to 
be Schur-concave, so they are entanglement monotones and may be used as alternative measures of 
entanglement. 
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I. INTRODUCTION 

Investigating properties of a quantum state it is useful 
to analyze its phase space representation. The Husimi 
distribution is often very convenient to work with: it ex- 
ists for any quantum state, is non-negative for any point 
a of the classical phase space O, and may be normalized 
as a probability distribution . The Husimi distribution 
can be defined as the expectation value of the analyzed 
state g with respect to the coherent state \a), localized 
at the corresponding point a fl. 

If the classical phase space is equivalent to the plane 
one uses the standard harmonic oscillator coherent states 
(CS), but in general one may apply the group-theoretical 
construction of Perelomov 0]. For instance, the group 
SU{2) leads to the spin coherent states parametrized by 
the sphere 5^ = <CP^ |3jB, while the simplest, degener- 
ated representation of SU (K) leads to the higher vector 
coherent states parametrized by points on the complex 
projective manifold <CP^-'^ |SlEl3- 

One may define 

the Husimi distribution of the A'^-dimensional pure state 
10) with respect to SU{K)-CS for any K < N, but it is 
important to note that if = any pure state is by 
definition S'C/(A^)-coherent. 

Localization properties of a state under consideration 
may be characterized by the Wehrl entropy, defined as 
the continuous Boltzmann-Gibbs entropy of the Husimi 
function ^J- The Wehrl entropy admits smallest values 
for coherent states, which are as localized in the phase 
space, as allowed by the Heisenberg uncertainty relation. 
Interestingly, this property conjectured first by Wehrl 8] 
for the harmonic oscillator coherent states was proved 
soon afterwords by Lieb 0], but the analogous result for 
the SU{2)-CS still waits for a rigorous proof. This un- 
proven property is known in the literature as the Lieb 
conjecture It is known that the SU {2)-CS provide 
local minimum of the Wehrl entropy [lol |. while a proof 
of the global minimum was given for pure states of di- 
mension iV = 3 and = 4 only 11, 12]. A generalized 
Lieb conjecture, concerning the Renyi -Wehrl entropies 
of integer order q > 2 occurred to be easier then the 



original statement and its proof was given in [Hill. It 
is also straightforward to formulate the Lieb conjecture 
for Wehrl entro py comp uted with respect to the SU{K)- 
coherent states but this conjecture seems not to 

be simpler than the original one. 

In this work we consider the Husimi function of an 
arbitrary mixed quantum state p of size A^ with respect to 
SU{N) coherent states, and calculate both its statistical 
moments and the Wehrl entropy. The difference between 
the latter quantity and the minimal entropy attained for 
pure states, can be considered as a measure of the degree 
of mixing. 

Analyzing pure states of a bipartite system it is helpful 
to define coherent states with respect to product groups. 
The Wehrl entropy of a given state with respect to the 
SU{N) X SU{N)~CS was first considered in the original 
paper of Wehrl Recently Sugita proposed to use 

moments of the Husimi distribution defined with respect 
to such product coherent states as a wa y to characterize 
the entanglement of the analyzed state [13 • 

In this paper we follow his idea and compute explic- 
itly the Wehrl entropy and the generalized Renyi- Wehrl 
entropies for any pure state \^) of the N x N composite 
system. The Husimi function is computed with respect 
to SU{N) X SU{N)-CS, so the Wehrl entropies achieve 
their minimum if and only if the state j^*) is a prod- 
uct state. Hence the entropy excess defined as the dif- 
ference with respect to the minimal value quantifies, to 
what extend the analyzed state is entangled. The Wehrl 
entropy excess is shown to be equal to the subentropy 
defined by Jozsa et al. |0| . Calculating the excess of the 
Renyi- Wehrl entropies we define the Renyi subentropy - 
a natural continuous generalization of the subentropy. 

Several different measures of quantum entanglement 
introduced in the literature (see e.g. 0, |23, 1^ and 
references therein) satisfy the list of axioms formulated 
in 19]. In particular, quantum entanglement cannot in- 
crease under the action of any local operations, and a 
measure fulfilling this property is called entanglement 
monotone. 

The non-local properties of a pure state of an A^ x A^ 
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system has to be characterized by iV" — 1 independent 
monotones . In this work we demonstrate that — 1 
statistical moments of properly defined Husimi functions 
naturally provide such a set of parameters, since they are 
Schur-concavc functions of the Schmidt coefficients and 
thus are non- increasing under local operations. 



momentum's z-component Jz with maximal eigenvalue, 
while the lowering operator reads Ji = J- — Jx — iJy 

To characterize quantitatively the localization of an an- 
alyzed state Q in the phase space we compute its Husimi 
distribution Hg with respect to SU{N)-CS and analyze 
the moments ruq of the distribution 



II. HUSIMI FUNCTION AND WEHRL 
ENTROPY 

Any density matrix g G (C^^^ can be represented by 
its Husimi function Hg{a), that is defined as its expec- 
tation value with respect to coherent states ja), 

Hg{a) = {a\Q\a) (1) 

In general one may define the Husimi distribution with 
respect to 5t/(iir)-coherent states with 2 < K < N , but 
most often one computes the Husimi distribution with 
respect to the SU{2) coherent states, also called spin co- 
herent states Ql^. Let us recall that any family of coher- 
ent states {|a)} has to satisfy the resolution of identity 

/ d^i{a)\a){a\ = 1 , (2) 
Jn 

where d^ia) is a uniform measure on the classical mani- 
fold Q.. In the case of the degenerated representation 
of SU{K) coherent states 5, 6, 7] this manifold, Hk — 
U{K)/[U{K - 1) X U{1)] = €P^-\ is just equivalent 
to the space of all pure states of size K. This complex 
projective space arises from the iiT-dimensional Hilbert 
space Tifc by taking into account normalized states, and 
identifying all elements of Tix , which differ by an overall 
phase only. In the simplest case of SU (2) coherent states 
it is just the well known Bloch sphere, ^2 = (DP^ = S'^. 

Analyzing a mixed state of dimensionality N we are go- 
ing to use SU (N) coherent states, which will be denoted 
by IttAf). With this convention the space is equivalent 
to the complex projective space (DP^^^, so every pure 
state is S'C/(iV)-coherent and their Husimi distributions 
have the same shape and differ only by the localization 
in n. 

The SU{N) coherent states may be defined according 
to the general group-theoretical approach by Perelomov, 
by a set of generators of SU{N) acting on the distin- 
guished reference state |0). We are going to work with 
the degenerate representation of SU (N) only, in which a 
coherent state may be parametrized by A'^ — 1 complex 
numbers 7i, 

|ajv) = |7i,...,7Ar-i):=e-^^-^^---e-^«--^"-|0) , (3) 

where operators J,; may be interpreted as lowering oper- 
ators . In language of the A^-level atom they couple 
the highest A^-the level with the i-th one (see e.g. 0). In 
the simplest case of SU (2) coherent states the reference 
state |0) is equal to the eigenstate \j,j) of the angular 



m,{g) = f dti{a) {Hg{a)f . (4) 

Here c?/i(a) denotes the unique, unitarily invariant mea- 
sure on (DP^-\ also called Fubini-Study measure. The 
measure d/i is normalized such that for any state g the 
first moment mi is equal to unity, so the non-negative 
Husimi distribution Hg may be regarded as a phase-space 
probability distribution. 

The definition of the moments vriq is not restricted to 
integer values of q. However, from a practical point of 
view it will be easier to perform the integration for integer 
values of q. But once ruq are known for all integer q, there 
is a unique analytic extension to complex (and therefore 
also real) q, as integers are dense at infinity. 
Another quantity of interest is the Wehrl entropy Sw^ 
defined as Q 

Sw{q) = - I dfi{a) Hg{a)\nHg{a) . (5) 
Jun 

Again, performing the integration might be difficult. 
But if all the moments niq are known in the vicin- 
ity of q — 1, one can derive Sw quite easily. Using 
dnydq = HnnH, one gets 

5^(,) = -lim^. (6) 
g^i oq 

The moments rUq of the Husimi function allow us to write 
the Renyi- Wehrl entropy 

Sw.qig) ■■= -^—\nmq{g) , (7) 

1 - q 

which tends to the Wehrl entropy Sw for g — > 1. As a 
Husimi function is related to a selected classical phase 
space, the Wehrl entropy is also called classical en- 
tropy 1^ Il6j |. in contrast to the von Neumann entropy 
Sn — —Trglng, that has no immediate relation to clas- 
sical mechanics. 

In section UTTI we consider how strongly a given state 
is mixed and in section IIVI we discuss the non-locality 
of bipartite pure states. For this purpose we pursue an 
approach inspired by the Lieb conjecture 1^, according 
to which the Wehrl entropy Sw of a quantum state is 
minimal if and only if g is coherent. In order to measure 
a degree of mixing of a monopartite mixed state of size N 
we therefore use the SU{N) coherent state, while for the 
study of entanglement of a bipartite pure N x N state 
we use the coherent states related to the product group 
SU{N) X SU{N). 
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III. MONOPARTITE SYSTEMS: MIXED 
STATES 

In this paragraph we wiU focus on mixed states g £ 
(^NxN acting on an iV-dimensional Hilbert space Hm- 
Such a state may appear as a reduced density matrix 
Qr, defined by the partial trace Pr(^) — Tr^|5')(^'| of a 
pure state l^*) G ®1-Ln of a bipartite system. One 
could also consider a system coupled to an environment, 
in which a mixed state is obtained by tracing over the 
environmental degrees of freedom. The von Neumann 
entropy Sn quantifies, on one hand, the degree of mixing 
of Qr, and on the other, the nonlocal properties of the 
bipartite state Y^). 

To obtain an alternative measure of mixing of g we 
are going to investigate its Husimi function defined with 

respect to the SU{N)-CS as ^(a) = {aN\Q\aN)- The 
moments niq of the Husimi distribution can be expressed 
as functions of the eigenvalues \i of g, which in the case 
of a reduced state of a bipartite system coincide with the 
Schmidt coefficients of the pure state \^). The deriva- 
tion of the explicit result in terms of the Euler Gamma 
function 



iv!r(g + i) 

T{q + N) 



with 



EA,, 
■ 



(8) 



N 



i=l 



n 



is provided in Appendix 1X1 Performing the limit (0 we 
find that the Wehrl entropy Sw equals the subentropy 
Q{g) up to an additive constant Cn, 

Swig) - Qie) + Cn . (9) 
The A^-dependent constant can be expressed as 



JV 



N 



■^{N + 1) - *(2) = 1/k 



(10) 



fe=2 



where the digamma function is defined by ^(x) = 
^''alz^^'* • '^^^ subentropy 

Q(g)--E =--QCx) (11) 

n (^^-^^■) 

was defined in an information theoretical context |l8l |. 
It is related to the von Neumann entropy S]\[{g) in the 
sense that both quantities give the lower and the upper 
bounds for the information which may be extracted from 
the state g The subentropy Q{g) takes its mini- 

mal value for pure states with only one non vanishing 



eigenvalue. Therefore we define the entropy excess AS' 
as 



AS{g) = Swig) ~ Swip^p) 



(12) 



where refers to an arbitrary pure state and the Wehrl 
entropy of an arbitrary iV-dimensional pure state is given 
by Swip-ip) — C'n- Hence the entropy excess, equal to the 
subentropy AS'(p) = Qig), is non negative and equal to 
zero only for pure states. 

As a byproduct we find a bound on Swio)- The suben- 
tropy Qig) is known not to be larger than the von Neu- 
mann entropy Sn E3- Using this fact and lO we find 
that Swig) ^ SNig) + Cn- As it is also known that the 
von Neumann entropy Sn is not larger than the Wehrl 
entropy Swig) 0j up with the following upper 

and lower bound on Swig) 



SNig) + Cn > Swig) > SNig) 



(13) 



IV. BIPARTITE SYSTEMS: PURE STATES 

Let us now focus on bipartite systems described by a 
Hilbert space Ti that can be decomposed into a tensor 
product H = Ha Hb of two sub-spaces. With a local 
unitary transformation Ui = hi a ®hiB any pure state 
l^*) e 7i can be transformed to its Schmidt form 



JV 



(14) 



where N =min(dim 7i^,dim TLb) and the real prefac- 
tors \i called Schmidt coefficients are the eigenvalues of 
the reduced density matrix gr- Thanks to the Schmidt 
decomposition we can assume dim T-La = dim Hb = N 
without loss of generality. Following an idea of Sugita 
|17| we consider the question whether the Wehrl entropy 
of a bipartite state can serve as a measure of entan- 
glement. He defined a coherent state {a^^'^'^) of an M- 
partite qubit system as the tensor product of M coher- 
ent states |a2^'') of single qubit systems. We general- 
ize this approach for bipartite systems, omitting the re- 
striction to qubits and calculate all moments and the 
Wehrl entropy of the respective Husimi functions. For 
a bipartite pure state \'^) G H, we use the tensor prod- 
uct of two SU(A^)-coherent states |a^^) = law)®^ S H 
to define a Husimi function. More explicitly one can 

write |q!]v ) = \ctN)A <^ \aN)B, with \aN}A G Ha and 
\o:n)b G Hb (we will drop the index referring to the 
subsystem, wherever there is no ambiguity). Such bi- 
partite coherent states were used already in the original 
paper of Wehrl . 

The Husimi function H oi a, bipartite state |^) is then 
given by 



i/*(aA, as) 



(^'KlaAr)^® |ajv)B)| 



(15) 
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By definition any product state is a SU{N) x SU{N) 
colierent state. According to the Lieb conjecture, orig- 
inally formulated for SU{2) coherent states, the Wehrl 
entropy is minimal for coherent states 9]. Hence it is 
natural to expect that the Wehrl entropy 



dfiA{a) / dfj.B{a) 
H-i,{aA, as) In iJ*(aA, as) 



(16) 



provides a measure of how "incoherent" a given state is. 
Due to the equivalence of "coherence" and separability 
we expect that Sw can also serve as an entanglement 
measure. In the following we discuss properties of Sw 
and show that it satisfies all requirements of entangle- 
ment monotones [T9ll20l|. 

There is a one-to-one correspondence between states 
and Husimi functions. Concerning entanglement, two dif- 
ferent states that are connected by a local unitary trans- 
formation are considered equivalent. As by construction 
the moments vriq of H are invariant under local unitary 
transformations, they reflect this equivalence and there- 
fore can be expected to be good quantities to character- 
ize the nonlocal properties of a bipartite state. For a 
pure state l^*) e T-La ^ Hb, there are — 1 indepen- 
dent moments, determining A^ — 1 independent Schmidt 
coefficients (one coefficient is determined by the normal- 
ization). In our case the moments read 



(17) 



where the integration is performed over the Cartesian 
product flpf X ftpf, being the space of all product pure 
states of the bipartite system. The proper normaliza- 
tion of dfiA and dfj,B assures that mi = 1. The re- 
quired A^ — 1 monotones can be provided by rUq with 
q = 2, . . . ,N. The moments can be expressed as a func- 
tion of the Schmidt coefRcients Xi, defined in eq (|14|l 



V r(g + A) 



(18) 



and are related to the monopartite moments (eq. ISJ by 
a multiplicative factor. The derivation of this result is 
provided in Appendix ^ Having closed expressions for 
the moments that can be extended to real g, one easily 
finds the bipartite Wehrl entropy 



The subentropy Q(A) takes its minimal value if and 
only if a given state is separable, i.e. if all but one 
Schmidt coefficients vanish. Therefore we define the bi- 
partite entropy excess AS'^(^) 



(21) 



that is vanishing if and only if \^) is separable. Fur- 
thermore, one can also consider the Renyi- Wehrl entropy 

In order to use the moments as entanglement mono- 
tones, it is conventional to rescale them such that they 
vanish for separable states and are positive for entangled 
ones 



M„ = 



1 



1 



-(Aiq.N - 1) 



(22) 



These quantities are analogous to the Havrda-Charvat 
entropy (also called Tsallis entropy) [23L |23| and in the 
limit g — > 1 one obtains the subentropy 



lim^ Af, = Q(A) 



RENYI SUBENTROPY 



(23) 



As discussed in previous sections, the excess of the 
Wehrl entropy A5 may be used as a measure of the de- 
gree of mixing for monopartie states, or degree of entan- 
glement for pure states of bipartite systems. Since the 
moments of the Husimi distribution are found, we may 
extend the above analyzis for the Renyi-Wehrl entropy 
Sw,q defined by Q. 

Considering, for instance, the case of mixed states of 
a monopartite system we use © to find the minimal 
Renyi-Wehrl entropy attained for pure states g^. 



1 ^^^A!r(g + i) 



(24) 



which for q^ \ reduces to (IIUII . In an analogy to © we 
define the Renyi-Wehrl entropy excess 



ASw,q{Q) — Sw,q{Q) — C^.q 



(25) 



Applying ||SJ| it is straightforward to obtain result in 
terms of the eigenvalues \i of the analyzed state q 



5'w('0) = - / d^i{aA) I rfM2(as) 
i/*(aA, as) In i/*(Q;A, as) • 



(19) 



Performing the limit ® one finds that Sw equals the 
corresponding monopartite quantity up to an additive 
constant 



Swi"^) = Q{\) + 2Cn ■ 



(20) 



/^Sw.qio) = 



1 



N 



A 



q+N-l 



N 



n (^«-^^) 



(26) 

This result shows that the excess of the Renyi Wehrl 
entropy Qq{g) — Qq{X) may be called Renyi subentropy 
since for q 1 it tends to the subentropy Hll|) . On 
one hand it may be treated as a function of an arbitrary 



5 




FIG. 1: Renyi entropy Sq (thin line) and Renyi subentropy Qq (thick hne) Une for the N = 2 probabiMty vectors as a function 
of the one independent variable x for g = l/2(a),g = l (b), g = 2 (c) and q = IQ (d). The maximum of the Renyi subentropy is 
increasing with q whereas the maximum of the Renyi entropy is independent of q. In the limit g — > oo both quantities converge 
to — InAmax, where Xmax is the largest component of the vector A. 
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FIG. 2: Renyi subentropy Qq as a function of two independent variables of the N = 3 probability distributions, q = 1/2 (a), 
q = 1 (b), g = 2 (c) and g = 5 (d). Gray scale is used to represent the values of subentropy: the higher q the larger the 
max;imum of Qq at the center A, = (1/3, 1/3, 1/3), which refers to max;imally entangled states. 
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FIG. 3: Rescaled moments Mq as a function of two independent variables for AT = 3 probability vectors, q = 1/2 (a), g = 1 
(b), g = 2 (c) and q — 5 (d). Labels at corners identify pure separable states. Note that in general Mq is not a monotonically 
increasing function of q. 



quantum state g, on the other it may be defined for an 
arbitrary classical probability vector A. 

In a sense Qq is a quantity analogous to the Renyi 
entropy Sq, which may be defined for an arbitrary prob- 
ability vector 



1 ^ 
1-0 



(27) 



In the case of quantum states Sq is defined as a function 



of the q-th moment of the respective state 

1 



1-q 



ln(Tr^«). 



(28) 



In the following we will sketch some properties of Qq{X) 
as a function of a classical propability vector. All the con- 
siderations arc also applicable in the quantum case, par- 
ticularly if we use Qq as an entanglement monotone using 
the moments discussed in the preceding section. Let us 
define two distinguished propability vectors which corre- 
spond to extreme cases: A* with X* = {i — 1, . . . , N) 
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describes the maximal random event, whereas with 
A^"'' = 6ij with j G [1, . . . , A^] describes an event with a 
certain result. For a bipartite quantum system a vector 
of Schmidt coefficients given by A* represents a maxi- 
mally entangled state, whereas As describes a separable 
state. 

The Renyi subentropy Qq{X) has the following proper- 
ties 

i. ) Qq{Xs) = for any q. 

ii. ) Qq takes its maximal value = 

r{q+i)m W for For q 1 one has 

Q^<^^ = Cn -\nN. 

iii. ) As tJq^'' = 1 one immediately has limq_^o Qq{^) = 

for any vector A. 

iv. ) For q — > 1 one obtains the regular subentropy, 

lim,_i Qq{X) = g(A). 

V.) For q ^ oo one gets lim^^oo Qq(A) = -InAmaa:, 
where Xmax is the largest entry of A. Hence this 
limit coincides with the limit q — s- oo of the Renyi 
entropy ^^(A). 

vi.) Qq is expansibile, i.e. it does not vary if 
the probability vector A is extended by zero, 
Qq{\l, . . . , Aat) ~ Qq{Xi, . . . , Aat, 0). 

A discrete interpolation between the subentropy Q and 
the Shannon entropy S was recently proposed in |25j |. 
Note that a continuous interpolation between these quan- 
tities may be obtained by the Renyi subentropy Qq for q 
increasing from unity to infinity combined with the Renyi 
entropy Sq for q decreasing from infinity to unity. It is 
well known that the Renyi entropy 5*^ is both a non- 
increasing function of q and convex with respect to q. 

We conjecture that the Renyi subentropy has opposite 
properties: it is non-decreasing as a function of q and it is 
concave with respect to q. See Fig. ^for some exemplary 
cases. Fig. ^ presents both quantities for N = 2, while 
Fig. 121 and Fig. |31 present the curves of isoentropy Qq 
and equal moments Mq obtained for = 3. 



or equivalently for the other considered quantities. The 
expression A ^ ^ means that A is majorized by ^, i.e. 
the components of both vectors listed in increasing order 
satisfy — J2i=i ^^r < j < N. In order to be 

Schur concave, ASw{X) has first to be invariant under 
exchange of any two arguments, which is obviously the 
case, and second it has to satisfy 



(Ai - A2) 



fdAS^ 



w 



dAS^ 



w 



V dXi 



dX. 



< 



(30) 



For convenience we will first consider Hq^N- Once we 
know about Schur convexity or Schur concavity of iiq.N 
one can easily deduce Schur concavity of all the other 
discussed quantities. The quantities ^J-q,N can be ex- 
pressed by the following representation (see eq. IBip 

f^q,N — ^(g"!^.^)'* /a dx{XxY ^ where dx is a short hand no- 
tation for dxi . . . dxN and the integration is performed 
over the simplex A containing all propability vectors x 
with -^i ~ 1- Making use of this relation and assuming 
Ai > A2 without loss of generality, we obtain 

dUq^N dfiq^N _ r(g + A^) 



dXi 



dXo 



r(q + i) 



dx iXxy-\xi-x2) . 

(31) 

The full measure of the propability simplex A can be 
divided into two symmetric parts, the first part specified 
by xi < X2, the second one by Xi > X2. After exchanging 
the labels xi and X2 in the second part, the first and the 
second part coincide and one gets 



dfiq^N d^q.N _ r(q + N) 



dXi 



8X2 r(q+l) 



dx{xi - X2){Xi - X2) 



I (Aixi -f X2X2 +y^XiXi] - {XiX2 + X2X1 +y^AtXj) j 

V i>2 i>2 I 



(32) 



Since x\ > X2 and Ai > A2, one has a::iAi+X2A2 > a;2Ai-|- 
X1A2. Therfore for q > 1 the integrand is non- negative 
and so is the integral. For < q < 1 the integrand is non- 
positive. Thus [iq,N is Schur convex for q > 1 and Schur 
concave for < q < 1. Using this, one easily concludes 
that 



(Ai -A2) 



d 
dXl 



d 
'3X2 



^^lq.N 



> 



(33) 



VI. SCHUR CONCAVITY 

So far we defined rescaled moments, Renyi entropy and 
Renyi subentropy. All these quantities are constructed in 
such a way that they are vanishing exactly for separable 
states. In order to legitimately use these quantities as en- 
tanglement monotones, we still have to show, that they 
are non increasing under local operations and classical 
communication, or analogously that they are Schur con- 
cave [20ll27l. i.e. 



X^^ 



ASwiX) > ASwiO 



(29) 



lim„ 



As ASw can be expressed as ASw — ^^^^^q^i ~dq~ 
Schur concavity of the subentropy can directly inferred 
from the corresponding properties of fiq^N- Also the 
rescaled moments (|22|) are Schur concave for all values 
of q, as 1 — q is negative for q > 1 where ^q^N is Schur 
convex. 

For the Renyi subentropy Qq with q 7^ 1 we have 



(Ai - A2) 



/ dQq dQq 



1 1 
^-qQq 



(Ai - A2) 



V^Ai 

/ d^J-q.N 

\ dX, 



8X2 

dy.q.N 
0X2 



(34) 
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FIG. 4: Renyi entropy Sq (dashed lines) and Renyi suben- 
tropy Qq (solid lines) as a function of q for two exemplary 
probability distributions. We have chosen A*' = 4 probability 
vectors with power law distributed components, pj oc j'^ with 
ft = 3/2 in a) and k = 3 in b). Stars at q = 1 represent Shan- 
non entropy (upper) and subentropy (lower). We conjecture 
that Qq is both non-decreasing and concave with respect to 
the Renyi parameter q. 



Qq is a positive quantity. Using Schur concavity of /iq.Ar 
for < g < 1 and Schur convexity for q > 1, we conclude 
that Qq is Schur concave for all positive values of q. 

Since we have shown that the subentropy, the rescaled 
moments and the Renyi subentropy vanish if and only 
if the considered state is separable and that these quan- 
tities are Schur concave, these quantities are entangle- 
ment monotones and may serve as legitimate measures 
of quantum entanglement. Let us emphasize that the 
monotones found in this work differ from the Renyi en- 
tropies |29j| and the elementary symmetric polynomials of 
the Schmidt coefficients [s^ and cannot be represented 
as a functions of one of these quantities. 



VII. OUTLOOK 

Defining the Husimi function of a given pure state of a 
bi-partite N x N system with respect to coherent states 
related to SU{N) x SU{N) product groups and comput- 
ing the Wehrl entropy allows us to establish a link be- 
tween the phase space approach to quantum mechanics 
and the theory of quantum entanglement. 

This approach may be considered as an example of a 
more general method of measuring the closeness of an 
analyzed state to a family of distinguished states. This 
idea, inspired by the work of Sugita |17(| , may be outlined 
in the following set-up: 

a) Consider a set of (pure or mixed) states you wish to 
analyze 

b) Select a family of distinguished (coherent) states |a), 
related to a given symmetry group and parametrized by 
a point on a certain manifold f2. This family of states 
has to satisfy the identity resolution, |Q;)(Q!|dri = I. 

c) Define the Husimi distribution with respect to the 'co- 
herent' states and find the Wehrl entropy for a coherent 
state, Srmn = S {\a)) 



d) Calculate the Wehrl entropy for the analyzed state, 

e) Compute the entropy excess AS = S{\^p)) — Smin, 
which characterizes quantitatively to what extend the an- 
alyzed state \tp) is not 'coherent'. 

For instance, analyzing the space of pure states of size 
N we may distinguish the spin SU{2) coherent states, 
or, in general the SU{K) coherent states (with K < N), 
parametrized by a point on CP^ and on CP^^^, respec- 
tively. In the case oi N x L composite quantum system 
we may distinguish the SU{N) (g) SU{L) coherent states. 
They form the set of product (separable) states and are 
labeled by a point on CP^^^ x CP^~^. For mixed states 
of size N we may select SU (N) coherent states, i.e. all 
pure states. In all these three cases we define the same 
quantity, entropy excess AS*, which has entirely differ- 
ent physical meaning. It quantifies the degree of non- 
5L/(iir)-coherence, the degree of entanglement and the 
degree of mixing, respectively, as listed in Tabled 

Describing the items a)-e) of the above procedure we 
have implicitly assumed that the Wehrl entropy is mini- 
mal if and only if the state is 'coherent'. This important 
point requires a comment, since the status of this as- 
sumption is different in the cases discussed. For the space 
of iV-dimensional pure states analyzed by Husimi func- 
tions computed with respect to SU{2) coherent states 
this statement became famous as the Lieb conjecture Q. 
Although it was proven in some special cases of low di- 
mensional systems 0, 0, , and is widely believed to 
be true for an arbitrary dimensions, this conjecture still 
awaits a formal proof. On the other hand, it is easy to 
see that the Wehrl entropy of a mixed state is minimal 
if and only if the state is pure, or the Wehrl entropy of 
a bipartite pure state is minimal if and only if the state 
is separable. This is due to the fact that in both cases 
the entropy excess is equal to the subentropy which is 
non-negative and is equal to zero only if the state is pure 
0. 

One of the main advantages of our approach is that it 
can easily be generalized to the problem of pure states of 
multipartite systems - see the last column of Table 1. It 
is clear that the entropy excess of j^*) is equal to zero if 
1 4*) is a product state, but the reverse statement requires 
a formal proof. Moreover, the Schmidt decomposition 
does not work for a three (or many)-partite case. Thus 
in order to compute explicitly the entropy excess in these 
cases, one has to consider by far more terms than in the 
bipartite case. In the special case of three qubits any 
pure state may be characterized by a set of five param- 
eters 31, 32J, and it would be interesting to express the 
entropy excess of an arbitrary three-qubit pure states as 
a function of these parameters. The second moment for 
this system was already calculated ^3 but expressions 
for other moments or for the Wehrl entropy are still miss- 
ing. 
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States 


Mixed 


Pure 


Pure 


Pure 


Pure 


Systems 


simple 


simple 


simple 


bipartite 


multipartite 


(Joherent 
states 


SU{N) CS 

|aiv> G CP^-' 


5(7(2) CS 

|a2) G CP' 


SU{K) CS 

Wk) g cp^-1 


SU(N) ® SU{N) CS 
|Qiv)®' GCP^'-' 


SU{N)®'^ CS 
|aiv)®^" G CP^"-i 


Minimal WeJirl 
entropy 


Smin = S{\aN)) 


Smin = 5'(|Q2)) 


Smin = ^dQif)) 


Smin = 25'(|a]v)) 


S'n.in = MS(|qjv)) 


Generic states 


p : TLn — > TIn 


\tP) G CP™-' 


IV-) G CP™"' 


|*> G (CP™-')''" 


1*> G (CP™-')''" 


WeJirl entropy 
excess AS 


S{p) - Cn 


Sim - C2 


^dV-)) - Ck 




^d*)) - MCn 


measures 
degree of 


mixing 


non SU{2) 
coherence 


non ^[/(i^) 
coherence 


bi-partite 
entanglement 


non- Al -partite 
separability 



TABLE I: Wehrl entropy excess in various set-ups. This quantity measures the 'distance' of a given mixed state from the 
subset of pure states; the 'distance' of a given pure state from the subset of all SU{K) coherent states; the distance of a given 
pure state of a composite system from the subset of all separable (product) states. 
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APPENDIX A: MOMENTS OF THE HUSIMI 
FUNCTION 

As already mentioned, every pure state belonging to an 
iV-dimensional Hilbcrt space TYat is a SU{N) coherent 
state and vice versa. Therefore, following we can 
parametrize all SU{N) coherent states by 



N-l 



1/2 

X, I |0) 



E 



XiC 



(Al) 



with < 



< 



1 - E = 



AT! 



N-1 

j=i+l ^3 



and d/i = 



— rii^i '^Vi- The states i = Q, . . . ,N — 1 
form an orthonormal basis, while |0) is the reference 
state. The Husimi function of a pure state Y^) 
with respect to SU{N)® SU{N)-coh.eT:ent states \a^^) = 
\oin)a ® |ajv)B reads 



N 

E 

1 



\/K\iA{aN\l^)A A{^J'\aN)A 
B{aN\v)B B{fJ-\aN)B , 



(A2) 



where is represented in its Schmidt basis, j'l') — 
X^iy V^W)a W)b- The q-th moment is than given by 



N 

E 



Ml ■ - ■Aig = l 



f^{JVi,...,JVN,Mi,...,MN) 



(A3) 



/<? = / dn{a){a\i^i) {m\a) {a\iy2) {fi2\a) ...{a\iyg) {fig\a) 



(A4) 

The integers Afi and Ai^ count how often the states {vi) 
and {fii\ appear in fq. Note that one has to integrate 
separately over both sub-systems. Due to the symmetry 
of the Schmidt decomposition with respect to the two 
subsystems, both integrations lead to the same result and 
it is just the square of the integral over fq that is entering 
Using the parametrization (|A1|) fq can be expressed 



niq 
as 



(Ni-Mi)ipi 



N-1 

^ ^^y~kT:^.-l^.+M^ (A5) 

where the integrations over performed in increasing 

order of i and the upper integration limit is given by 
^(max) _ ^ _ ^^-^^ Xj. It is useful to perform the (pi 
integrations first. They lead to S_\f.^Mi terms, so that one 
gets 

fq = ^' ( n "'^^^ ^f'^M.M.j 



N-l 



(A6) 



In order to perform the integrations over Xi we need to 
define an auxiliary function 



9q 



(/3,A^) = f n / dx. 1 (/5-E -.)^-^"- 

(A7) 

with x["^'^^\p) = P^J2f=i+i ^j- the following we will 
show that for any integer q 



9 



holds. According to the definition (|A7|) . gq(f3, N) satisfies 
the following relation 



To end, we will show by induction that for any positive 
integer q the quantity fiq,N can be written as 



Jo 

(A9) 

Making use of 



dxi xf' {P-xif-^^ 



(9 + 1)! ' 
we immediately have (|A8p for N ^ 2 

rl3 



(AlO) 



(All) 



V- Af-2 



A/'i=0A^2=0 A/'n-i=0 

(A17) 

As a starting point for this line of reasoning we are going 
to demonstrate that H-i^n = 0. For N = 2 this equality 
can be checked by direct calculation. For N > 2 we have 



A 



N-3 



A,: 



N 

^N-2 



+ 



(A18) 



Assuming that (jA8|) holds true for gq{(3, N —1) and mak- 
ing use of (jA9|) one gets 



gq{P,N) ^ Jl^ dxN^i x^"^''{f3 - XN~i) 

(g+Ar-2-A^„-l)! 

Using further l|A10|l one has 



\q+N-2-J\rN 



(A12) 



a(3N) - 3,+N-l ^fN-ll{q + N-2-^^M-ly. 

3^^^''^> - P {q + N-l)\ 



{q + N -2-Nn-i)\ 
Finally one ends up with 



gq{P,N)=l3'^+''- 



-i (nt7'M!)(g-Efj'M)! 

{q + N-l)\ 



(A14) 

i.e. (|A8p also holds for gq{(3,N). Finally for /? = 1, we 
get 



{q + N~iy. 



(A15) 

Now we have = I],yi...iy,=i '^i/i ■■•'^i',/? ' where the 
(^n^.m^-terms assure that there are only integer powers 
of the Ai/;. Since fq depends only on the integer pow- 
ers Afi of the Xi,. , we need to count how many terms 
with fixed powers occur. Using simple combinatorics one 

can see that there are (^ (jj"-i(A/;)!)(g_^"-iA/-^)! ^ terms 

containing the expression A^^A^^ . . . A^"^^A^ ^i=i 
Thus finally we get 

= EE- E 



A"i=0A/'2=0 A/'iv-i=0 



^2 ... '^N^i 



(A16) 



Since 



= 1 



j^Tj^, we can make use of 

the following relation U.=^,.X^,i^^->^^) = 

fJ-q.N-iiXi, ■ ■ ■ , Xj^i, Xj+i, . . . , Xn), in which the right 
hand side vanishes by assumption. Therefore we end up 
with 



N 



M-i,Af - Xj ^ 



A, 



N-3 



~i nf=i,fc5^i(-^» - Xk) 



(A19) 



We still have the freedom to choose the index j. As the 
result does not depend on the choice of j, both sides of 
the equality have to be zero. Now we can come back to 
the proof of (|A17p . For TV = 2 it is just a straight forward 
calculation to show that 



E^i^r^ 



Afi Xl+' 



Ai — A2 A2 — Ai 



(A20) 



which can be checked by multiplying both sides with Ai 
A2. Assuming that HA17() is true for — 1 we get 



rUq^N 



(A21) 



Making use of the assumption (|A17|I one obtains 



JV-l 

E 



xr' 

A? 



E 



1 rijlij/iiAj - Xj) 



N 



\N-2 \q+l \q+l 

rrn,tM^.(A,-A,) A.-A^ 

Now one just has to make use of /i-i,Ar = and one im- 
mediately gets that the assumption ljA17|) holds true for 
N . This result together with (jA16|l completes the proof 
of Eq. ijSJl for integer positive integers q. Since integers 
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are dense at oo, we conclude that there is a unique gen- 
erahsation to real q. Thus eq ||SJ) is also valid for real 
q. Eq (|18|l for the moments of the Husimi function of a 
monopartite system differs by a proportionality constant 
only and its proof is analogous. 



APPENDIX B: INTEGRAL REPRESENTATION 
OF THE MOMENTS 



In order to prove that 

T{q + N) 



T{q+1) 



(Bl) 



where the integration is performed over the propability 
simplex A, we consider the function hq^Nia, b) of two real 
variables a and b defined as 

"-9,^^ = "FT — TTT / dxn-i ... dxi 
r(g + 1) Jo Jq 

n— 1 n— 1 

+ (1 - a - ^ x,)K + bf . (B2) 



Now we are going to show by induction that h can be 
expressed as 



^ ((l-a)A. + b)^+"-^ 



(B3) 



where we were using the assumtion for hq,N-i{ci + 
x„_i, 6 + a;„_i A„_i. Performing the integration one gets 



hq^N — 



T{q + N) ((l-a)A,+5) 



q+n—l 



r(g + i) 



E 



n,^.(A.-A,) 
^ (A^-i-A.)n,^,:(A.-A,: 



We are done, if we manage to show that 



N-2 

E 



1 



(A^-i-A.)n;;7(A.-A,) 
1 



riili^iAAf-i - Aj) 



or eqivalently 



P 



N~2N-2 . . 

\ ^ 1 r Ajv_i — \j 



n 



Xi — A," 



1 



(B5) 



(B6) 



(B7) 



This is a polynomial of {N — 2)-nd order in Aat-i. If we 
find N — 1 different values for Ajv-i such that p equals 
1, the polynomial p has to be identically equal to unity. 
Inserting A^ (fc = 1, . . . , iV — 1) one gets 



For 2 it is straight forward to check that the as- 

sumption EH holds. For > 2 we obtain 



hq^N — 



nq+N) 

T{q 



N-2 



1) h n,#. 



1 



(Ai — Xj) Jo 



dxN-l 



((1 - a - a::Ar_i)AAr_i + b + xn-iXn-i) 



q+N~2 



,(B4) 



Xh — A, 



(B8) 



Thus p = 1, which completes the proof of eq. (|B3|) . Set- 
ting a = 6 = 1 in eq. ljB2p and eq. ljB3|l , we immeadiately 
get eq. (|BT|| . 
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